Controlling soliton excitations in Heisenberg spin chain through magic angle 
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We study the nonlinear dynamics of collective excitation in a iV-site XXZ quantum spin chain, 
which is manipulated by an oblique magn etic field. We show that, when the tilted field is applied 
along the magic angle 60 = ± arccos \/l/3, the anisotropic Heisenberg spin chain becomes isotropic 
and thus an free propagating spin wave is stimulated. And in the regime of the tilted angle larger 
and smaller then the magic angle, two types of nonlinear excitations appear, which are bright soliton 
and dark soliton. 
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I. INTRODUCTION 

A crucial challenge for quantum information process- 
ing (QIP) is perfectly transmitting a quantum state from 
one place to another as well as storing the information 
of a quantum state. Recently, as much attentions have 
been paid to the solid-based quantum computing, there 
are many proposals for quantum state transfer and stor- 
age, using quantum spin systems as quantum data bus 
and quantum memory p| H, H, [1] • 

In other hand the Heisenberg chain is a typical spin sys- 
tem in condensed-matter physics. It has attracted con- 
siderable attentions for a long time in describing various 
magnetic properties of materials. As a strong correlated 
system [E |(| 0, HI, @, it can display rich nonlinearities 
induced by inter-spin interaction. The typical nonlinear 
phenomenon is the solitary wave, that is a wave packet 
propagate without either energy loss and wave packet 
spreading. In this sense it is quite natural to consider 
the role of soliton wave in quantum state transfer and 
quantum information storage. 

Generally there are two major approaches to theoret- 
ically study the solitary excitation in the one-dimension 
(ID) Heisenberg chainjlll, El El EI El EE El El El • 
1) each spin is represented as a Bloch vector, and a soli- 
tary excitation is found by investiga ting the nonlinear dy- 
namics of this classical variable [12|, LL3|- 2) Another pro- 
cedure is to employ the boson mappings of spin operators 
via the Holstein-Primakoff transformation and Jordan- 
Schwinger realization By making use of the spin- 
coherent-state representation, a solvable nonlinear differ- 
ential equation can be deduced from the Heisenberg equa- 
tion for the mapping boson operators, therefore a soli tary 
excitation is implied in the system (HI fTH El. fT7l ITl. fl9|. 

Our present investigation is motivated by the well 
known knowledge in nuclear magnetic resonance (NMR) 
that the effective spin-spin interaction can be controlled 
by the angle of the external applied magnetic field [Iol|. 
Here, to some extent, the inter-spin coupling can be can- 



celed effectively by the magic angle . Nowadays together 
with the spin echo technology, magic angle manipulation 
has became a necessary tool for controlling the interact- 
ing spin system to be free of decoherence due to the inter- 
spin couplings. In this paper, we apply a magnetic field 
on a N-site XXZ spin chain, which is rotated around the 
y-axis by an amount of magnitude 6, and study how to 
manipulate a soliton excitation by this oblique field. It 
is known that the ferromagnetic order is completely de- 
termined by the direction of the external magnetic field. 
Indeed, our investigation shows that the switch between 
the bright soliton and the dark soliton is controlled by 9, 
whose tangent is defined as the ratio between the z com- 
ponent and the x component of the magnetic field; and at 
the turning magnitude referred to as the magic angle Oq, 
the original anistropic chain becomes isotropic, therefore 
only an ideal spin wave is stimulated as an effective plane 
wave. Here, we employ the second approach mentioned 
above, but go beyond the spin wave approach by consid- 
ering the nonlinear effect of the collective spin excitation. 
Generally, the spin wave approach is regarded as a mean 
field method with a given order parameter. Due to the 
ignorance of the nonlinear fluctuation, the spin wave ap- 
proach does not cover the quantum fluctuation in the 
nonlinear regime. However our approach overcomes this 
disadvantage. 
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The remainder of this paper is organized as follows. 
In section mi we present our model - a Heisenberg spin 
chain with anisotropic coupling in an oblique magnetic 
field. The approach we use is in sections IHII and IIVI 
Here the quasi-classical equation of motion is obtained, 
which is a nonlinear Schrodinger (NLS) equation. In sec- 
tionfVl the time evolution of nonlinear excitations is stud- 
ied. The bright soliton and dark soliton can be excited 
as 8 shifted to the left or the right of the magic angle. In 
section IVI1 we give a remark about the relationship be- 
tween the soliton wave propagation and perfect quantum 
state transfer. 
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II. A SPIN CHAIN MODEL 

We consider a one-dimensional XXZ spin chain with 
N spins in an homogeneous external magnetic field 

B = B (e x sin 6 + e z cos 8) (1) 

Denote the spin operator at the jth site by Pauli spin 
operator S*, Sj, and 5J. Then Hamiltonian H of this 
system reads 

N 

H = - J E (Sj^i+i + s - s :!- ' + A ^; +1 ) + B s ( 2 ) 

i=i 

where S = • Sj is the total spin operator. Three types 
of interactions are included in Eq.@: 1) the isotropic 
part of nearest neighbor exchange interaction J x = J y = 
J > 0. 2) the anisotropic part of nearest neighbor ex- 
change interaction, which is characterized by a dimen- 
sionless parameter A = J z /J. 3) the external homo- 
geneous field B, which is specified by its components 
B x = £?sin# and B z — BcosO. Here the magnitude 
of the magnetic field B = WB^. + Bl , and the oblique 
angle 9 = arctan(_B z / B x ). 

Obviously, the last term is proportional to the total 
spin z-component Sjfotal = Ei^?; which is conserved 
and thus has common eigenstates with the total Hamil- 
tonian. Correspondingly, the Hilbert space of the system 
can be decomposed into a direct sum of numerous sub- 
spaces V(M) specified by the total spin number M along 
z-axis. It had been proved [2(|[2l[ that, on a finite sim- 
ple cubic lattice, the ground state of the XXZ model is 
nondegenerate in an subspace V(M). In particular, its 
global ground state ^>o( A) is just the ground state of 
the model in the subspace V(M = 0) When the lattice is 
finite, the ground state energy Eq(A) and the spin cor- 
relation function (S?S?) are analytical with respect to 
A. 

Here, we assume that the external magnetic field is 
much larger than the inter-spin interaction. Therefore 
we switch to a new "reference frame" with new Z— axis 
along the quantized direction of spin along the effective 
field B. The corresponding — > Lj transformation for 
each spin at site j reads 

L* = Sf cos 0-S? sin^, 

q = s?, (3) 

L) = £| cos 6 + sine. 

Therefore Hamiltonian in Eq. @ is decomposed into a 
direct sum of the SO(3) irreducible tensors with respect 
to Lj , i.e, 

M=2 

H = BL Z + H m (4) 

M=-2 

which is written according to the irreducible representa- 
tion DM (L = 2, 1, 0)of SO(3) group. Here, L z = £\ L) 



is the third component of the total angular momentum. 
With 5= (A -1), 

j 

- 6 - J (3 cos 2 Q-1)Y,L]L) +1 (5) 

3 

belongs to the zero rank representation D>- ' of 50(3), 
while 

H ±1 = |jsin2^ (L' s Lf +1 +LfL' j+1 ) (6) 

3 

and 

H±2 = -~J sin 2 9 Y. L t L t+i ( 7 ) 

3 

belong to the 1st and 2nd rank representation £>W and 
I)' 2 ' respectively. The Mth order tensors Hm satisfy 

[L Z ,H M ] = MH M , [L z , H ] = 0. (8) 

From the point of view of perturbation theory, the 
eigenvalues of the Hamiltonian are those of BL Z at the 
Oth order, i.e. L z = M, which means the ground state 
has all spins parallel to the field and it is not degenerate. 
However, the excited state, where one or two spins are 
flipped with respect to the ground state, are degenerate. 
Therefore, when the perturbation theory is applied, we 
diagonalize Hamiltonian H in each subspace of a given 
M . which means that at lowest order, Hamiltonian Hq is 
kept. Thus, the Hq plays an indispensable role to govern 
the dynamics of the spin chain. It can be found directly 
from Hamiltonian Hq that the spin-spin interaction is 
controlled by the direction of th e applied magnetic field. 
When 6 = 9q = ± arccos y/l/3, the original anisotropic 
chain becomes isotropic. 

We further explain the above argument from the point 
view of the representation theory of SO (3) together with 
rotating wave approximation. Sandwiched by the com- 
mon states \ J,M') and |J, M") of Hq and L z , the equa- 
tions (JSj) lead to 

{J,M'\[L Z ,H M ]\J,M") 
= (M 1 - M"){J,M'\H M \J,M") (9) 
= M(J,M'\H M \J,M") 

here M — M' — M" for the nonvanishing matrix element 
(J, M'\Hm\J, M"). In the interaction picture, the off- 
diagonal elements ( J, M'\Hm\J, M") are fast changing 
with high frequency, therefore the terms with \M\ = 1 
and \M\ = 2 in the total Hamiltonian can be ignored [lol ]. 

III. QUASI-CLASSICAL MOTION EQUATION 
FOR FERROMAGNETIC SPIN CHAIN 

With the above considerations, the Hm (M = ±1, ±2) 
are the first and second order tensor operators that trans- 
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form according to the representation D^ =1 ' 2 ^ of 5*0(3) 
group. In large external field limit, Hamiltonian Q is 
reduced to 



mation should be taken 



H = BL 7 



Ho 



(10) 



In this section, we begin with the above Hamiltonian to 
discuss the soliton excitations. 

To describe the spin excitation, we introduce the boson 
excitation by the Holstein-Primakoff transformation [13] , 



2S — OjCLj, 



at a j 



s, 



(11a) 

(lib) 
(11c) 



where the annihilation operators dj and the creation 
operators a\ satisfy the boson commutation relation 



di,a 



char- 



= Sij. The number operator n, 

acterizes the spin deviation from its maximum value S of 
if . Usually, for a Heisenberg ferromagnetic system with 
inter-spin couplings — JLj-L i+ i (J > 0), the spontaneous 
symmetry breaking will happen along the direction of the 
external magnetic field (says along z-axis), which is sub- 
sequently assumed to approach zero. The ground state 
implies a ferromagnetic order with all spins along the +z 
or — z direction. Therefore the magnetization comes into 
being, which is defined as the non-vanishing average of 
L z = Yli^i- The governing equation for the nonlinear 
excitations on ground states can be obtained by the ex- 
panding V25 - hi ~ V2S(1 - hi /{AS)) to a series over 
the low excitation {hi) -C 2S. We first expand (|lla|) and 
(lllb[) . and keep the terms in the first order of (hi) 



V2Sa\ ( 1 - 



(12a) 
(12b) 



Consequently, beyond the spin wave approximation, 
the low energy effective Hamiltonian (TlT)|) is achieved as 



H = —cqS (^jO'j+i + h-cj + B n j 

3 3 

+ir (''<"/"•• 1 + n A«]+i + h - c ) 

3 

- (co + ci) ^rij (rij+i - 2S) , (13) 



where rij = dlcij, cq — J (l + 5 sin 2 9/2) and c\ = 
5 J (3 cos 2 6>- 1) /2. 

As we should emphasize that, the external magnetic 
field has been assumed B > in the above discussions. 
If B < 0, another form of Holstein-Primakoff transfor- 



L) 



a\\ 2S — a\a 



3 

= s 



3 3 ' 



a]aj, 



(14a) 

(14b) 
(14c) 



but the main result of this paper is same. So without loss 
of generality, B > is assumed in the rest of this paper. 

Hamiltonian (fT3"|) characterizes the low energy nonlin- 
ear property of ferromagnetic spin chain in an oblique 
magnetic field. The corresponding Heisenberg equation 

ifr^- = — c S (&j+i + a,j-i) + (2S (cq + ci) + B) a,j 



Co 

' 4 

- (co + ci) d 3 n i±1 



(15) 



contains various nonlinear couplings. They lead to differ- 
ent nonlinear "phases" , which are represented by various 
types of soliton excitations. 



IV. CONTINUUM FIELD APPROACH FOR 
NONLINEAR COLLECTIVE EXCITATION 

In this section, we make use of the continuum field the- 
ory to study the nonlinear excitations in the ferromag- 
netic spin chain. To this end we consider the spin wave 
approach with some nonlinear corrections. First let us 
introduce the p— representation (or call the Glauber co- 
herent state representation [23[) defined by the product 
of the multi-mode coherent states \a) — J\i \ a i), where 
each component \a/) is the eigenstate of the annihila- 
tion operator dj, i.e., hi \a/) = on \ot\), and on is the 
coherent amplitude. Since coherent states are normal- 
ized and overcompleted, the field operator sandwiched 
by | a) can be represented only with their diagonal ele- 
ments. Thus, we only need to consider the diagonal part 
of Eq. (|15[) . which are enough to describe the nonlinear 
dynamics without any help of off-diagonal elements. The 
p-representation of non-linear equations (| 15[) reads 

iK-^- = -cqS (ctj±i - 2atj) + (2Sc x + B) otj 

+ -J ( 2 Kl 2 a j±i + K±i| 2 «j±i + a* ±1 a jaj 



(Co + Ci) OSj \aj±l\* 



(16) 



From then on, the spin dynamics is expressed in terms 
of the c— number equation in Glauber's coherent-state 
representation. However Eq. (fll))) is difficult to solve due 
to its nonlinearity and discreteness. Since we are look- 
ing for excitations with a length scale much larger than 
the lattice constant, we take the long- wave approxima- 
tion, that is, the continuum field theory approach will 
be employed in the large N limit. By assuming that the 
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coherent amplitude is continuum in space, the discrete 
variables oti{t) can be replaced by a mean field cp(z, t), 
i.e. at (t) — > tp (z,t). Correspondingly, the difference be- 
come the following differential 



Q!i±l 



az 



1&_ 

2~d7 2 



(17) 



where the lattice constant is assumed to be unity. By 
keeping the derivation terms to the second order d 2 /dz 2 , 
the system of equations (fl6l) become a field equation with 
self-coupling terms 

d d 2 
ih—ip + c S-^(p + 2ci(p\if\ 2 = Vip, (18) 

which is a typical NLS equation in a constant potential 
V = 2S Cl + B. 

To separate the rapid oscillation in the inhomogeneous 
nonlinear equation (|18[) . we rewrite ip(z,t) as 



(19) 



where £ = Xy/K/JfioS), \ = (2Sci + B) /h, and (/>(£, t) 
is the slow varying envelope. Then the standard form of 
the NLS equation is obtained 



*0t + 0« + 2^- 10| 2 . 



0, 



(20) 



which has the soliton solution. Here, the strength of 
nonlinear term 



Cl 



-J (3 cos 2 1 



1 



(21) 



is determined by spin-spin interaction parameter J and 
the angle 9 of the obliqu e ma gnetic field. We notice that 
when 9 = 6o = arccos y/l/3, the nonlinear term in the 
equation of the motion disappears, therefore Eci. (f20]) be- 
comes a standard wave equation, which gives a complete 
set of plane waves solution. In this case, the collective ex- 
citations are spin waves. 9q is named as the magic angle. 
It shows that the magic angle changes the anisotropic 
system to the isotropic one. 



V. NONLINEAR TIME EVOLUTION OF 
LOCALIZED MAGNETIZATION 

For this ferromagnetic spin chain, the interaction pa- 
rameter J is positive. Thus the sign of coefficient c\ is 
determined by the angle of the oblique magnetic field 9. 
And nonlinear property of this system is strongly related 
with the sign of the nonlinear term. With these different 
physically accessible parameters, the nonlinear equation 
(|2U)1 can possess the bright-soliton and the dark-soliton 
solution. Next, we use the inverse scattering method 
[2~il |25| to obtain the different solitons of this system. 

When ci > 0, that is S (3 cos 2 9-1) > 0, the single 
bright soliton solution of Eq. (TIB]) is obtained as 



ip (x, t) = Ae^-^sech 



Cl 



Atl— l — (x — xo — vt) 
c b 



, (22) 



with v — viy/ cqS/H, 7 = Hv/2cqS and 



to 



Cl 



(2S - A 2 



B 
J 



hv 2 

4c^S 



(23) 



where v\ is an integral constant. The coefficient v — 
vi y/coS/H describes the velocity of bright soliton travel- 
ing to the right. The positive coefficient A characterizes 
the size of the bright soliton. The coefficient Xq denotes 
the initial position of the bright soliton. Parameters vi, 
A and Xq are determined by the initial state. The bright 
soliton solution (f2"2"|) depicts a wave packet traveling in 
the continuous background. 

When ci < 0, that is S (3 cos 2 9—1) < 0, we can get a 
single dark-soliton solution of Eq. (fT8|) 



Lp (x, t) = A'e 1 ^'*-"'*) tanh \a\ (x - x' - v't) 



hv 



with v' = v'yy/coS/h, i - 2cqS 
. r.<7- 

n 



and 



4c 5 



(24) 
(25) 



where v[ is an integral constant, v' = v^y/coS/H charac- 
terizes the velocity of dark soliton traveling to the right. 
The positive coefficient A' characterizes the size of the 
dark soliton. The coefficient x' Q denotes the initial posi- 
tion of the dark soliton. Parameters v', A' and Xq are 
decided by the initial state. The dark soliton solution 
describes a localized dip in the continuous background. 
Fig. [1] numerically illustrates the bright and the dark soli- 
ton obtained from Eq. (f2"2"]) and respectively in com- 
parison with the results with many extra nonlinear terms 
in the following discussion. 




50 100 

X 

FIG. 1: Bright soltion (a) obtained from Eq. J22J) for 9 = 0.1. 
Dark soltion (b) obtained from Eq. flMJ) for 9 = 1.5. Other 
parameters are taken as follows S = 10, B = 100, J — 1,5 = 
0.1, x = x' = 0, A = A' = 1, vi = v[ = 5 and t = 3. 

Furthermore, to test the rationality for neglecting the 
terms H±i and H±2, we deal with Hamiltonian ([4]) by 
the approach presented above. The Holstein-Primakoff 
transformation (jTTJ) is first employed to write Hamilto- 
nian (|U) in terms of bosonic operators 
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H = -coS (a]a j+1 + 6jaj +1 J + B ^ nj - (c + a) ^ n s (n 3+1 - 25) 

3 3 3 



2c 3 5 ^ ( -— ^otnj-at +1 + a]a] +l n ]+1 ^j + a)a] +1 + h.c.) , (26) 



where c = J (l + 5 sin 2 0/2), c x = J5 (3 cos 2 9 - l) /2, magnetic spin chain in an oblique magnetic field. The 
c 2 = J<5sin26>/4 and c 3 = J 5 sin 2 9/4. Hamiltonian ([26]) corresponding Heisenberg equation 
characterizes the low energy nonlinear property of ferro- 



daj 
ih—r- 
dt 



-c Sdj±i + (B + 25 (c + ci)) a,j — (c + Ci) ajUj ±1 + (2rijdj±i + TijiiOjii + a] ±1 ajOj 
-C2V2S ^j±i + + %±iSj — 25 + Tij + -ajSj^ 

-c 3 25 (a] ±1 ~ Jg n 3^±i ~ ^ d ]±i n J±i ~ ^%%%±i) 



r 



(27) 



contains rich nonlinear couplings and thus can predict 
various nonlinear phenomena. 

Further we use the continuum field theory approach, 
that is, expressing the spin dynamics according to its cor- 
responding continuous c-number equation in Glauber's 
coherent-state representation. Then a similar field equa- 
tion is obtained as 



ih^-ip + c S-^tp + 2cnp \ip\ 



dt 



dx 2 



{2Sc 1 +B)i P + 2c 2 v / 2S | - \tp\ 



7<P ~S 



-C3 



(*Stp* - 3<p* 



(28) 



Obviously, the effects on H±i and H±2 is kept in the 
above derivation. Eq. ([2"8")) is similar to the NLS equa- 
tion, but more nonlinear terms are involved, therefore 
the analytic solution can not be got. 

To display the effects of extra nonlinear terms, we 
numerically investigate the time evolution of dark and 
bright solitons under the action of differential equa- 
tion dHJ) in Fig. [2] and Fig. [3 Figures [2] and [3] show 
that a bright and a dark solition can be excited when 9 
takes value on the left and the right side of the magic an- 
gle respectively. Comparing the numerical results shown 
in Fig. [2] and [3[ it is found that the bright soliton is more 
stable than the dark soliton solution. This is because the 
long time evolution of the dark soliton in Fig. [3] become 
a profile with spatial oscillation. While the bright soli- 



ton will almost keep their shapes, the dark soliton will 
disappear after a long enough evolution. Thus the bright 
soliton is more easily detectable in practice. From Fig. [21 
we found that the closer to the magic angle the tilted 
angle is, the wider the soliton becomes. 

In Fig. [4] we plot the time evolution for the stable 
bright soliton with respect to different values of the di- 
mensionless parameter A = B/J. It can be found that 
the extra nonlinear terms has no effect on solitons when 
A = 100. Such numerical analysis confirms that the 
soliton-like wave can be excited under an appropriate 
regime from A = 10 to 1000 approximately. Therefore 
in this regime, it is reasonable to discard the terms H±i 
and H±2 in Eq. Q. When parameters A goes to infi- 
nite (see the line of A = 5000 in FigfJJ, the interaction 
term Hq can be completely ignored, then Hamiltonian 
of the system approximately described by BL Z . In this 
situation, the system is completely polarized along the 
direction of the magnetic field and form a background 
with symmetry breaking which is displayed by the dot 
line in Fig. [4j The numerical result that a very sharp 
peak is localized around x = justifies our observation 
from the physical intuition. 

There are some experimental data on quasi-one- 
dimensional ferromagnetic chains [13, [H, HtJ [3(| l3l| . 
which may be used to test our predictions in an in- 
direct way. For the spin chain made of the material 
CagCoa-^Fe^Oe (x=0, 0.1, 0.2) (Ref. [13), the exchange 
constant J decreases from 4.39 to 8.13K and the mag- 
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FIG. 2: Numerical time evolution of a bright soliton under 
the act of Eq. flH]| for 9 = 0.1 (a) and 9 = 0.9 (b). Other 
parameters are taken as follows S — 10, B = 100, J = 1 and 

8 = 0.1. It indicates that a bright soliton is excited when 

9 < O . 




50 100 



FIG. 3: Numerical time evolution of a dark soliton by the 
influence of extra nonlinear terms in Eq. (|28[) for 9 = 1.0 
(a) and 6—1.5 (b). Other parameters are taken as follows 
S = 10, B = 100, J = 1 and 8 = 0.1. It indicates that a dark 
soliton is excited when 9 takes value on the right side of the 
magic angle. 



netic field strength B < 10T (~ 33.6if). The solitary- 
excitation is possible since the dimcnsionless parameter 
A = BJ J 10. However, for the three ferromagnetic 
chains [H, [2^, [3(| [Mj], which made of the following ma- 
terial: (1) [(CH3) 3 NH]FeCl3-2H 2 (Ref. [H,[2|), where 
J ~ 17 AK,B < 16T(~ 21.5K); (2) CoCl 2 -2D 2 0(Ref. 
[13), where J ~ 2A75K and B < 6T (~ 12. IK); 
(3) CoCl 2 -H 2 0(Ref. HI), J - 18.3K and B < 6T 
(~ 12.1K), the bright soliton may be excited, but noises 
exist as shown in Fig. U with the dimensionless parame- 
ter A = B/J w 1. 



VI. 



REMARKS AND CONCLUSION 



Before concluding this paper, let us discuss the phys- 
ical meaning of the solitary wave of magnetic excitation 
in a N-site XX Z quantum spin chain from the point 
view of quantum information processing. We first note 
that the multi-mode coherent states \a) — Y[j \<Xj) repre- 
sent an inhomogeneous collective excitations distributed 
around the spin chain. This collective excitation has a 
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FIG. 4: (Color online) the time evolution of a bright soliton 
under the effluence of extra nonlinear terms in Eq. (|28|l with 
respect to different values of the dimensionless parameter A = 
B/J at 9 = 0.1, S = 10,8 = 0.1, t = 3. 
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spin representation by 



1 



= Q L 1 



(29) 



with Q = 1/y/S -Lf + l and 10)^ = \S, -S) j being the 
lowest eigenstate of the on-site spin Iy. We also notice 
that the coherent state denotes a superposition of various 
spin states \S, m s ) -, i.e., 



^C„(a 4 )(gL+)"|5,- 



|-B n (aj) 15,71- SV 



(30) ' 

where C„(ai) = exp(— |ay| /2)a!y 3 /(nj-!) and B n (oti) can 
be obtained from C n (ai) directly. During the propaga- 
tion of soliton, the narrow traveling wave packet does 
not spread, therefore it behaves like a "flying qudit" (a 
d = 2S level system). 

When S = 1/2, the qubit is localized at the jth site 
with a superposition state 



1 1 

2' 2 



(31) 



which can be used to encode quantum information as 
usual. During the propagation, the wave function nearly 
keeps its spatial shape all the time. From an mathemat- 
ical point of view, the spatially non-spreading properties 
of the carrying excitation wave is very crucial for quan- 
tum state transfer from one location to another with high 



fidelity. It seems the bosonic excitations obey the bosonic 
commutation relations only in the large S limit, but the 
spin wave approach can still work well for S = 1/2 in the 
condensed matter physics. In this sense we can suppose 
our above arguments available. 

In summary, we have studied solitary magnetic excita- 
tion in a iV-site XXZ quantum spin chain as well as 
how to use an oblique magnetic field to create differ- 
ent types of solitons. Through a mean field approxima- 
tion beyond the usual spin wave approach, we obtain the 
quasi-classical motion equations for nonlinear evolution 
of the Heisenberg spin system. We show that the switch 
between the bright and the dark soliton is controlled by 
the angle of the magnetic field, whose tangent is defined 
by the ratio between the z component and the x com- 
ponent of the magnetic field. And at the magic angle, 
the system go to the isotropic Heisenberg model, hence, 
only an ideal spin wave is stimulated. We also remark 
the possibility for solitons to play the role of "fly ing qu- 
dit" based on the well-known results in Ref. [26fl that a 
no-spreading wave-packet behaves like a "flying qubit" . 
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